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I. INTRODUCTION 
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The chiral phase transition at finite temperature T and/or chemical potential fj, is studied using 
the QCD-like theory with a variational approach. The "QCD-like theory" means the improved 
ladder approximation with an infrared cutoff in terms of a modified running coupling. The form 
of Cornwall-Jackiw-Tomboulis effective potential is modified by the use of the Schwinger-Dyson 
equation for generally nonzero current quark mass. We then calculate the effective potential at finite 
T and/or /i and investigate the phase structure in the chiral limit. We have a second-order phase 
transition at Tc — 129 MeV for fi = and a first-order one at /ic ~ 422 MeV for T = 0. A tricritical 

- - - point in the T-^ plane is found at T = 107 MeV, /i = 210 MeV. The position is close to that of the 

^5 random matrix model and some version of the Nambu-Jona-Lasinio model. 
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I At zero temperature and zero (baryon number) density, the chiral symmetry in quantum chromodynamics (QCD) 
. is dynamically broken. It is generally believed that at sufficiently high temperature and/or density the QCD vacuum 
^ ■ undergoes a phase transition into a chirally symmetric phase. This chiral phase transition plays an important role 
OO in the physics of neutron stars and the early universe and it may be realized in heavy-ion collisions. At finite 
temperature the lattice simulation is powerful to study the chiral phase transition at finite temperature {T ^ 0). It 
, is now developing also for finite chemical potential (/z ^ 0). However, efi^ective theories of QCD are still useful for 
' various nonperturbative phenomena including the phase transition. 

' Recently, it has been argued that importance of a study of the phase structure, especially a position of the tricritical 
point, has been pointed out in Ref. [|l|. The Nambu-Jona-Lasinio (NJL) model |^ in which the interaction is induced 
by instantons and the random matrix model have shown almost the same results concerning the tricritical point. 
'"q^: It is also interesting to study the possibility of color superconducting phase at high baryon density [2,4-8]. However 
I we may neglect this phase in the high temperature region where we found the tricritical point. In this paper, we 
Qh! concentrate on the chiral phase transition between SU{Nf)L x SU{Nf)]i and SU{Nf)L+R using the effective potential 
^ ' and the QCD-like theory. One usually studies the phase structure of QCD in terms of the Schwinger-Dyson equation 
(SDE) or the effective potential [9-13]. However, the use of the SDE only is not sufficient for its study in particular 
when there is a first order phase transition; then, we use the effective potential. The QCD-like theory provided with 
the effective potential for composite operators and the renormalization group is successful to study the chiral symmetry 
breaking in QCD [|l^,^. This type of theory is occasionally called QCD in the improved ladder approximation. The 
cd phase diagram in the QCD-like theory has been studied in Refs. p[ pT|Jl^ . However, the position of the tricritical 
point is largely different from that obtained from the NJL model and the random matrix model. 

In this paper we use a modified form of the Cornwall-Jackiw-Tomboulis (CJT) effective potential [|l^ which is 
convenient for a variational approach. The formulation is given for the case where the chiral symmetry is explicitly 
broken at zero temperature and density. We, then, consider the CJT effective potential in the improved ladder 
approximation at finite temperature and/or density. Being motivated by Refs. [§,p|, we re-examine the chiral phase 
transition and phase structure in the chiral limit. 

This paper is organized as follows. In Sec. II we formulate the effective potential for composite operators and 
extend it to finite temperature and density. In Sec. Ill we first determine the value of Aqcd by a condition = 93 
MeV at T = /i = and then calculate the effective potential at finite T and/or /i numerically. Using those results, we 
study the phase structure in the T-fj, plain. Sec. IV is devoted to conclusion. We fix the mass scale by the condition 
Aqcd = 1, except for Sec. III. 
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II. EFFECTIVE POTENTIAL FOR QUARK PROPAGATOR 



A. CJT effective potential at zero temperature and density 



At zero temperature and zero density, the CJT effective potential for QCD in the improved ladder approximation 
is expressed as a functional of S{p) the quark full propagator 



V[S] = Vi[S] + V2[Sl 

Vi[S] = J ^ Tr [lniS,\p)Sip)) - S^\p)Sip) + l] , 
d'^p d'^q 



V2[S- 



-C2 



{2tj:)H (27r)4i 



: f{p, g)Tr {l^S{p)^,S{q)) D^^^ [p ~ q), 



(1) 
(2) 

(3) 



where C2 = {N^ — l)/(2iVc) is the quadratic Casimir operator for color SU{Nc) group, So{p) is the bare quark 
propagator, g^{p, q) is the running coupling of one-loop order, D^^^{p) is the gluon propagator (which is diagonal in 
the color space) and "Tr" refers to Dirac, flavor and color matrices. The two-loop potential V2 is given by the vacuum 
graph of the fermion one-loop diagram with one gluon exchange (see Fig. 1). 

After Wick rotation, we use the following approximation according to Higashijima and Miransky 



g'^iPE, qs) = 0{pE - qE)9^iPE) + 0{qE - PE)9^{qE)- 



(4) 



In this approximation and in the Landau gauge, no renormalization of the quark wave function is required [ p9[ and 
the CJT effective potential is expressed in terms of S(p_e) the dynamical mass function of quark: 



v[npE)] 

+4 

y2[S(p£)] =-6^2 



Vi[nPE)]+V2[nPE)l 

-2 



(27r)4 m2(A)+p| 
^ d*PE S(pg)(I](pg)-m(A)) 
(27r)4 J:Hpe)+P% 
^ '■^ d'^PE d^qs g^iPE,qE) 



Y,{pi 



(2^)4(27r)4 {pE 

nqE) 



qs) 



■i:Hpe)+pI mqE) + qV 



(5) 



(6) 



(7) 



Here, an overall factor (the number of light quarks times the number of colors) is omitted and to(A) is the bare quark 
mass. In the above equations we temporary introduced the ultraviolet cutoff A in order to make the bare quark mass 
well-defined. 

The extremum condition for V with respect to S(p_e) leads to the following SDE for the quark self-energy 



1:{pe) = m(A) + 3C2 



(i^qE 9^iPE,qE) ^(qE) 



(27r)4 (pE-qE)^ J:^qE)+ql- 



(8) 



In Higashijima-Miransky approximation, since the argument of the running coupling has no angle dependence, we first 
perform the angle integration. As a result, we understand that the procedure is achieved equivalently by replacing 
iPE - Qe)'"^ by 9{pe — qE)ip%)~^ + 0{qE — PE)iqE)~^ Elq. (H). Then we can reduce Eq. (||) to the following 
differential equation 

(4^)2 



^^(Pe) 

and the two boundary conditions 



■Pi 



1 d^pE) 
3C2 p%dpl {A{pe) dpi 

1 dS(pB) 



^(pe) 



Mpe) dpi 

V{pe) d^{pE) 
^{pe) dpi 



= 0, 



PE=^ 



PE=A 



m{A), 



(9) 

(10) 
(11) 
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where the functions 



Pe 

and 

d 



T^KPe) = 2 — (12) 



are introduced. 

Substituting Eqs. (|) and @ into Eqs. (|) and (0), we obtain 

'•^ d^PE , ^Hpe)+pI 



nnPE)] = -2j i^m 



m2(A) +p| 
^ dV S(pg)(S(pg)-m(A)) 
(27r)4 mPE)+pl 

^\pe)+pI 
(27r)4 ''m2(A)+p|; 



[E(pb) - m(A)] 



3C2 J (2^)4 ^ ^-^-^ ^ p'^dp'^y^iPE) dpi 

= 1 ^^iPE)+pi 

J (27r)4 m2(A)+p| 
2 /"^' . . I f d 



where we used a partial integration in the last line and 

Vs ^ FiA) ~ F{0), 

F{pe) ^ ^ nPE) - m(A)] (15) 
30*2 A(p£;) dpi; 

Hereafter we consider the effective potential in the continuum limit (A oo). Let us begin by evaluating F(A) 
using the running coupling 

-1, N 27r2a 24 

9{Pe) = -, 2 ' "^TiT? ^ — ' 16 

mpj; llA'c — 2n/ 

and the corresponding asymptotic form of the mass function 

/I 9 \ —a/2 

^lnp^\ ^^,,2, a/2-1 



E(p^) ^ ™(A) j +_(lnpi)'^/-^ (17) 

where n/ is the number of flavors which controls the running coupling. Throughout this paper, we put Nc — Uf ~ i, 
namely a = 8/9. The parameter a is related to the order parameter of the chiral symmetry {qq) as 

a = (18) 

When the chiral symmetry is exact, i.e., m(A) = 0, using Eqs. (|l^ ) and (|l7|), we can easily show that F(A) vanishes 
in the continuum limit, i.e., limA^oo F{^) = 0. As for F(Q), since we introduce infrared finite running coupling and 
mass function in Eqs. (|20| ) and (^), we can set F{0) = 0. After all, in the continuum limit, we get T/g = and the 
modified version of the CJT effective potential is obtained as IlBJ2fl] 



d^PE , ^^{pe)+pI 



pI 
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We can also show that Vs = 0, namely Eq. (|l^) holds for nonzero bare quark mass 
A few comments are in order. 



(1) The extremum condition for Eq. (19) with respect to T,{pe) leads to Eq. (|9|) which is equivalent to the original 
equation (||) in Higashijima-Miransky approximation apart from the two boundary conditions. We will take account 
of these conditions when we introduce the trial mass function. 

(2) In the chiral limit, Eq. (^9|) is the same as the expression given in Refs. ||l^,|2^. However, even if the chiral 
symmetry is explicitly broken, we can use the same expression for V pl|] . We do not require the finite renormalization 



adopted in Ref. [15 



Now we are in a position to introduce a modified running coupling and a trial mass function. We use the following 
QCD-like running coupling |14| 



ln(p| + p\) ' 



fiPE) = ^Z7:^r-Za\^ (20) 



where is a parameter to regularize the divergence of the QCD running coupling at p = 1(Aqco). This running 
coupling approximately develops according to the QCD renormalization group equation of one loop order, while it 
smoothly approaches a constant as p\ decreases. 

Hereafter we consider the chiral limit; i.e., the m(A) = case. Corresponding to the QCD-like running coupling, 
the SDE with the two boundary conditions suggests the following trial mass function [Oj 



2 , 2 a/2-1 

Pe ~^ Pr 

where a is the same as before 



S(p^) = ^— ^ [HpI + pD] "'^-^ , (21) 
Pe ~^Pr 



Using Eqs. (20) and (pi[), we can express y[I](p£;)] as a function of a the order parameter. A further discussion 
of the CJT effective potential and the dynamical chiral symmetry breaking in QCD-like theory at zero temperature 
and density can be found in Refs. 0,^. 



B. Effective potential at finite temperature and density 

In this subsection we discuss the effective potential at finite temperature and density. In order to calculate the 
effective potential at finite temperature and density, we apply the imaginary time formalism 

, oo 

-^fiVi) -^T /(a;„ + , {n e Z), (22) 

n— — oo 

where ujn — (27^ -I- 1)t:T is the fermion Matsubara frequency and /i represents the quark chemical potential. In 
addition, we need to define the running coupling and the (trial) mass function at finite T and/or /i. We adopt the 
following real functions for T>T,^iji) and T,t,^{p) 



27r^a 1 

•Ct.mIW = 1 . 2 I 2T 2' (23) 



St,.(p) = ,,2^" , ,2 H-I+P'+Pir"^' ■ (24) 
P ^ Pr 



In Eq. ( p3| ) we do not introduce the /i dependence in 'Dt,^i{p)- The gluon momentum squared is the most natural 
argument of the running coupling at zero temperature and density, in the light of the chiral Ward-Takahashi identity 
p3| . Then it is reasonable to assume that T>t.^{p) does not depend on the quark chemical potential. 

As concerns the mass function, we use the same function as Eq. (^ ) except that we replace p^ with a;„. As already 
noted in Sec. II A, the quark wave function does not suffer the renormalization in the Landau gauge for T — fi = 0, 
while, the same does not hold for finite T and/or /i. However, we assume that the wave function renormalization is 
not required even at finite T and/or /i, for simplicity. 

Furthermore, we neglect the T-fj, dependent terms in the quark and gluon propagators which arise from the per- 
turbative expansion. We expect that the phase structure is not so affected by these approximations. 

Using Eqs. ( p3| ) and (p^, it is easy to write down the effective potential at finite temperature and chemical potential 
(see Appendix). Assuming the mean-field expansion, the effective potential can be expanded as a power series in a 
with finite coefficients a2n {T, fi) 
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V{a;T,fi)= a2 (T, fi)a' + ai{T, fi)a^ + ■ ■ ■ . (25) 

Once we know the value of (Tmin the location of the minimum of V, we can determine the value of (qq) using the 
following relation 



{qq) = -TJ2 I J^'^^StAp)^ (26) 



(27r)3 



where ST,fj.{p) is the quark propagator at finite T and/or fi in our approximations and "Tr" refers to Dirac and color 
matrices. However, in this paper, we still determine the (qq) through the relation {qq) — — (3/27r^a)f7mm. We have 
confirmed that this relation works well even at finite T and/or ^. 

III. CHIRAL PHASE TRANSITION AT HIGH TEMPERATURE AND DENSITY 

In our numerical calculation, as mentioned before, we put Nc = nj = 3. Furthermore, since it was known that 
the quantities such as {qq) and /^r are quite stable under the change of the infrared regularization parameter , 
we fix tR = In(p|j/Ap(,j3) to 0.1 and determine the value of Aom by the condition /tt = 93 MeV at T = /_* = 0. We 
approximately reproduce /jr using the Pagels-Stoker formula [|24| : 

and obtain Aqcd = 738 MeV. The value of Aqcd is almost the same as the one obtained in the previous paper [|l2[ in 
which we used Eqs. (|) and (|). 

A. T / 0, = case 

Fig. 2 shows the T-dependence of the effective potential at /i = 0. We can realize that (Jmin the minimum of the 
effective potential continuously goes to zero as temperature grows. Thus we have a second-order phase transition at 
Tc = 129 MeV. Fig. 3 shows the temperature dependence of —{qq)^^^- 

B. T = 7^ case 

Fig. 4 shows the /x-dependence of the effective potential at T = 0. For small values of ^, the absolute minimum 
is nontrivial. However we find that the trivial and the nontrivial minima coexist at /i = 422 MeV. For larger values 
of /i, the energetically favored minimum move to the origin. Thus we have a first-order phase transition at /ic = 422 
MeV. Fig. 5 shows the chemical potential dependence of —{qq)^^^- The chiral condensate vanishes discontinuously at 
M = Mc- 

C. T / 0,/i / case 

In the same way as the previous two cases, we determine the critical line on the T-fi plane (see Fig. 6). The position 
of the tricritical point "P" is determined by the condition 

a2{Tp,^ip) = ai{Tp,^p) ^Q, (28) 

in Eq. (p5|). Solving this equation, we have 

{Tp,^p) = (107,210) MeV. 

We have varied tn from 0.1 to 0.3 in order to examine the tp dependence of the position of P. As a result, for instance, 
we have 

{Tp,fip) = (104, 207) MeV for tp = 0.2, 
= (101, 208) MeV for tp = 0.3. 

We note that the value of Aqcd has been determined at T — fi = hy the condition — 93 MeV for each value of 
tp. Thus we confirmed that the position of P is stable under the change of tp. 
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IV. CONCLUSION 



In this paper we studied the chiral phase transition at high temperature and/or density in the QCD-hke theory. 

We extended the effective potential to finite T and fi and studied the phase structure. We found the second-order 
phase transition at Tc = 129 MeV along the = line and the first-order phase transition at /Zc = 422 MeV along 
the T = line. We also studied the phase diagram and found a tricritical point P at (Tp, fip) = (107,210) MeV. 
Phase diagrams with similar structure have been obtained in other QCD-like theories pJlSl . As concerns the position 
of the tricritical point, however, our result is not close to theirs. Let us consider the reason why our model gives the 
different result. In Ref. they used the momentum independent coupling and the mass function without logarithmic 
behavior. The values of Tc and /ic of Ref. [|l3| are about the same as ours. However, the position of the tricritical 
point is in the region of small fj,. The discrepancy may arise from the fact that: (1) They did not use the variational 
method, but numerically solved the SDE; (2) The treatment of the gluon propagator at finite T and/or /i is different 
from ours. Our result is rather consistent with that of the NJL model [g[ and the random matrix model They 
obtained 



Tp - 100 MeV , 3fip - (600 - 700) MeV. 

Recently it was pointed out that the values of T and fj, accomplished in high-energy heavy-ion collisions may be close 
to the tricritical point and it may be possible to observe some signals ||l|. Thus it is significant that three different 
models show almost the same results. 

Finally, some comments are in order. In this paper, we modified the form of the CJT effective potential at 
T = fi = using the two representation of the SDE. Our formulation of the effective potential is entirely based on 
the Higashijima-Miransky approximation. It was known that the approximation breaks the chiral Ward-Takahashi 
identity. Therefore, it is preferable to formulate the effective potential without this approximation. However it seems 



that the results do not depend on choice of the argument momentum |23 2^. Moreover, the treatment of the quark 
and the gluon propagators at finite T and/or /i is somewhat oversimplified in the present work. We would like 
to consider the wave function renormalization and more appropriate functional form for I?T,/i(p) and Tit,^{p)- By 
including a finite quark mass, we can study a more realistic situation where the chiral symmetry is explicitly broken. 
In the studies of SDE with the finite quark mass, it was known that there is a difhculty in removing a perturbative 
contribution from the quark condensate p9| , p6[ |. In the effective potential approach, however, we are free from such 
a difficulty. The study of the phase structure with the finite quark mass is now in progress [ pl| . Furthermore, we 
also plan to study the quark pairing including a color superconductivity [2,4-8] and a "color-flavor locking" |Q (for 
Nc — Nf — 3 case), in the QCD-like theory. 



APPENDIX: 



In this appendix, we show the effective potential explicitly. In the first place, we consider the case of zero temperature 
and finite chemical potential. 
Using Eq. (p^), we obtain 



d*PE ^^ ^t.m(p) + (P4+^/^)'+p' 



1 

47r3 



(2^)4 
In 



(P4 + 2^)2 -I- p2 



where the imaginary part of Vi is odd function of p^; therefore it has been removed from Eq. ( |A1[ ) and 



(Al) 



dpi 



d\p\ p2 



(A2) 



In Eq. (19), we carry out the momentum-differentiation and, then, use Eqs. ( |23| ) and (p4[). V2 is obtained as 

y _ 4a^ [ {pI+p^?Mp1^p^+pI)Y -^ 



37r3C2a 7p {pI + p^ + p^) ln(p| + + p|) + pj + 



1 



(p2+p2+p|)3 



HpI + p^ +pI) + i-- 



(A3) 
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At finite temperature and chemical potential, the integration in Eqs. (Al) and (A3) is replaced by the sum over 
the Matsubara frequencies. 
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FIG. 1. Two-particle irreducible graph which contributes to V2. 
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FIG. 2. The effective potential at finite temperature and zero chemical potential. V is defined hy V = 24-!T^V and all 
quantities are taken to be dimensionless. The curves show the cases T/Aqcd = 0, 0.15, 0.2. 
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FIG. 3. The temperature dependence of —{qqY^^ at /^ = 0. 
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FIG. 4. The effective potential at finite chemical potential and zero temperature. The curves show the cases ///Aqcd = 0, 
0.5, 0.6. 
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FIG. 5. The chemical potential dependence of —{qq)^^^ at T = 0. 
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FIG. 6. The phase diagram in the T-fi plane. Solid line indicates the phase transition of second-order and points indicate 
that of first-order. The point P is the tricritical point. 
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